TUM-HEP-603/05 



Leptonic Dark Energy and Baryogenesis 
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We consider a baryogenesis scenario, where the difference of baryon (B) and lepton (L) number is 
conserved in such a way that the B—L asymmetry in the standard model sector is compensated by 
an asymmetry of opposite sign stored in the dark energy sector. Therefore, we introduce a toy-model 
in which a complex quintessence field carries a B—L asymmetry at late times. We determine the 
produced baryon asymmetry in the visible sector for a large range of initial conditions and find it 
easy to achieve a value of the observed order of magnitude. While the size of the produced baryon 
asymmetry depends on details of the underlying infiationary model, it turns out to be independent 
of the reheating temperature in many cases. We also discuss possible sources of instability like the 
formation of Q-Balls. 

PACS numbers: ll.30.Fs, 98.80.Cq 



I. INTRODUCTION 

One of the most intriguing questions in cosmology and 
particle physics concerns the origin of the observed asym- 
metry between baryons and anti-baryons [1]. If we as- 
sume that our universe underwent an inflationary phase, 
any initial asymmetry has been diluted during inflation 
and some baryogenesis process must have generated the 
observed asymmetry afterwards. Frequently it is stated 
that any baryogenesis process has to fulfill the three 
Sakharov conditions: a) B-violation, b) C- and CP- 
violation, and c) departure from thermal equilibrium [2] . 
Since only B—L is conserved in the standard model (SM), 
most baryogenesis models introduce new B and/or L vi- 
olating interactions at some high energy scale. Typical 
scenarios of this kind include heavy particle decays or 
make use of scalar field dynamics. An example for the 
latter is the Affleck-Dine mechanism [3], where complex 
scalar superpartners of fermionic fields acquire baryonic 
and/or leptonic charges due to B—L violating interac- 
tions in their potentials. These asymmetries are later on 
transferred to the SM fermions. However, it is also pos- 
sible to create a baryon asymmetry in the SM particle 
sector within a B-conserving theory. All one has to do 
is to store a B (or B—L) asymmetry of opposite sign in 
some hidden sector, see e.g. Ref. [4]. The case that this 
sector is dark matter and that it is completely made up 
of antibaryons has first been been considered in Ref. [5] , 
where the ratio of the energy density of SM particles and 
dark matter leads to a prediction for the masses of the 
constituents of dark matter. Another example for a B—L 
conserving baryogenesis scenario is found in Ref. [6] . 

Recently, another topic in cosmology has also drawn 
a lot of attention since one has found evidence [1, 7-10] 
that the expansion of the universe is accelerating. One 
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origin of this acceleration might be some unknown en- 
ergy form called dark energy (DE), see Refs. [11-14] for 
some recent reviews. Quintessence, a homogeneous scalar 
field rolling down a suitable chosen potential is a popular 
candidate for DE [15, 16]. In most quintessence models 
the scalar field is real-valued, however, the variant of a 
complex scalar has also been considered [17, 18]. The 
major difference results from the possibility to assign a 
global conserved U(l) quantum number to the complex 
scalar field, i.e. it can carry a charge. Here, a possible 
connection to baryogenesis emerges. 

Since scalar fields could be responsible for the observed 
baryon asymmetry as well as the acceleration of the uni- 
verse, it seems natural to consider the possibility that 
a scalar field could be connected to both phenomena. 
While Ref. [18] mentions several possible realizations for 
such an idea, including the one of hiding anti-baryons in 
the vacuum in a baryon-symmetric baryogenesis model, 
this paper introduces a first model (to our knowledge) 
that draws a direct connection between the charge of a 
complex field responsible for dynamical DE and the ob- 
served baryon asymmetry in our part of the universe. For 
various other possible connections of baryogenesis and a 
DE scalar field or the vacuum see also Refs. [19-25]. 

In the presented scenario B—L is conserved and a cor- 
responding asymmetry is stored in the DE sector com- 
pensating for the observed baryon asymmetry and possi- 
ble other asymmetries in hidden sectors. Therefore we in- 
troduce a simple toy-model, in which an additional scalar 
field mediates between a complex quintessence field (car- 
rying lepton number) and the fermionic sector. At early 
times relative phase differences in the initial conditions 
will lead to lepton asymmetries in both condensates. The 
additional field eventually decays and transfers its asym- 
metry to the fermionic sector, where sphalerons [26, 27] 
create a baryon asymmetry. The new field also acquires 
a huge mass due to the large vacuum expectation value 
(VEV) of the quintessence field and hereby effectively 
decouples the fermionic and the DE sector. 

The paper is organized as follows: in Sec. II we define 
the Lagrangian of the model and discuss the initial con- 
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ditions and parameters. Wc also explain the evolution 
of tiie scalar fields and the order of events that lead to 
the final baryon asymmetry. A quantitative treatment 
of the scenario can be found in Sec. Ill, where we give 
analytical estimates for the final baryon asymmetry of 
the scenario while we also present numerical results for 
a large range of initial conditions. In Sec. IV we discuss 
several possibilities for the needed coupling between the 
SM and the newly introduced scalar sector. In Sec. V 
we study the stability of our model concerning particle 
processes and inhomogeneities in the scalar condensates, 
that might lead to Q-Ball formation. Finally, we present 
our conclusions in Sec. VI. 



II. THE SCENARIO 

Our scenario begins at the end of inflation, which 
ensures a flat and homogeneous Friedmann-Robertson- 
Walker space-time. Ignoring SM particles and right- 
handed neutrinos for the moment the Lagrangian L of 

the system is given by 

+ \{d^^r{d>^x)^\^l\x? (1) 

- iAi|,^nx|2-iA2(</.V' + h.c.), 

with the quintessence field (/> and a mediating field x- 
For stability reasons we require Ai > A2 > 0. While the 
X field will later on be coupled to the fermionic sector, 
which results in an additional term in its equation of mo- 
tion, we take the behavior of the quintessence field to 
be completely determined by the classical equations of 
motion derived from this Lagrangian. This is equivalent 
to the assumption that all interactions of arc already 
included in its effective potential V(<^) and the given cou- 
plings to the X field. This will be further motivated in 
section V. 

One can see that the system has a global C^(l)- 
symmetry. Later on, we will choose the fermionic cou- 
plings in such a way that we can identify the correspond- 
ing conserved quantity with lepton number L and as- 
sign both fields the charge L = —2. Therefore, the pos- 
tulated L''(l)-symmetry for the complete model is just 
B—L, which is already inherent in the SM. The lepton 
number density stored in the quintessence condensate is 
in this case given by 

n^ = -i{r<l>-(t>*4>) = -2\ct>\^d, (2) 

where 9 is the phase oi (p = jc/ile'^. Then the lepton 
number density per comoving volume has the form 

A, = -2\^f() (3) 

where a{t) is the cosmic scale factor and hq = a(to)- 
Furthermore, the quantities and A^ corresponding to 



the mediating field x = |x|e"^ with the phase function <t 
are defined analogously. 

For the quintessence potential we choose a negative 
exponential form [28, 29], 

V{\ct>\) = Vo (e-^^l-^l/^P' + fce-^^l-^l/^P'), (4) 

where Mpi = I/^/SttG is the reduced Planck mass. This 
is a typical quintessence potential which leads to a light 
quintessence mass and extremely small higher deriva- 
tives at late times. Also, for typical parameter values 
(^1 = 0(10), ^2 = 0{1), < 1) such a potential yields a 
tracker behavior, where the energy density of the quint- 
essence field follows the background energy density. The 
dark energy and cold dark matter densities arc directly 
related by a constant ratio of order S/^f [29] during mat- 
ter domination in these models. The second term in the 
potential yields the transition to the late-time acceler- 
ation. Since our work concerns only early times, we 
set /c = for the rest of this paper. Also many other 
models lead to equivalent potentials at early times, see 
e.g. Refs. [30, 31). Typical values of the quintessence vac- 
uum expectation value at late times are of the order of 
the Planck scale and beyond. 

Our scenario starts at the end of inflation with an os- 
cillating inflaton and Hubble rate iJinf- The amplitude 
of the inflaton oscillations will first decay due to expan- 
sion of the universe and later on due to the decay of the 
inflaton condensate during reheating {H = H^). Obvi- 
ously these conditions are fulfilled by most inflationary 
scenarios and a dependence of our scenario on the partic- 
ular inflationary model only comes into play when quan- 
titatively calculating the produced amount of the final 
baryon asymmetry. 

At the beginning {t = to = 0) both flelds have initial 
VEVs of respective absolute amounts (po, Xo and a rela- 
tive phase ao- Similar initial conditions were considered 
in Refs. [32, 33] (see also Refs. [34, 35]). We assume the 
VEVs to be of such kind that the main part of the energy 
density is still given by the inflaton. 

The equations of motion resulting from the La- 
grangian (1) are given by 

4> + 3H^ = ~2^-X,\xU-X2rx\ (5) 

X + 3Hx + rx = -(M' + Ai|(/.nx-A2X>', (6) 

where we added an additional damping fac- 
tor r = g^m^l {Stt) due to the fermionic couplings 
of X [36, 37], with g-^ being the sum of the squares 
of the coupling constants of these interactions and 
= (yu^ -|- Ai|0p)-^/^ being the effective mass of x-^ In 



Even though different directions of x might have different masses, 
the definition of F is perfectly fine, as long as wc arc only con- 
cerned with orders of magnitude and Ai and A2 are not fine-tuned 
(see also Sec. III). For the same reason possible decays of x mto 
</) can also be neglected (see Sec. V). 
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order for the interaction of 4> and x to have a noticeable 
effect we require /i^ ^ Ai|0p at all times. We also 
stress, that whenever this damping term comes into 
effect the relation |^/0| <C (t> is fulfilled and we can 
consider as approximately constant. 

Since the last term in brackets of Eq. (f) can be 
written as ^A2|(/)p|x|^ cos(2a), with the phase differ- 
ence a = 6 — a, we see that both fields will develop 
non-zero phase-velocities if the initial conditions ful- 
fill sin(2ao) ^ 0, corresponding to a spontaneous CP- 
violation, see also Ref. [32]. 

In the case of negligible fcrmionic couplings (F = 0) 
lepton number is conserved in the scalar field sector and 
as long as the homogeneous condensates do not decay 
the phase space rotation of (j) and x will always fulfill the 
constraint Aff, = —A^. 

For sizeable fermionic couplings of x this simple re- 
lation no longer remains valid. Using Eq. (3) and the 
equations of motion (5), and (6) we now find 

^(A^ + A^) = -r^^, (7) 

which shows that the explicit development of the fields 
has to be considered, which in return depends on their 
initial conditions at the end of inflation. For the moment 
we will fix them in a way where one can easily sec how a 
baryon-asymmetry in the fermion sector can arise. How- 
ever, the system is unlikely to be restricted to these initial 
values and we will consider a much broader parameter 
range in the following section. With this in mind, we put 
A2 of order one tenth, all other coupling constants of or- 
der one and both initial field values of the order of Hi^i at 
the end of inflation. Additionally, we choose Vb in such a 
way that i?iifM|i > F(0 w Hi^f) > H^^^fMpi, which en- 
sures that the main part of the background energy den- 
sity at the beginning of the scenario {pq"^ = SHf^fMp^) 
is given by the infiaton, while the behavior of cf) is still 
governed by V{(j)) and not by its mixing terms with x- 

With the above mentioned, </> will now rapidly increase 
right from the start for typical parameter values in the 
quintessence potential. Since the mixing terms arc of 
the order of the Hubble scale, the phase velocities of 
both fields also become sizeable. When grows big 
enough such that F w J? the x-condensate starts to decay 
due to its fcrmionic couplings and its oscillation ampli- 
tudes rapidly decrease. Since these interactions are lep- 
ton number conserving the decay also mediates a lepton- 
number transfer to the fcrmionic sector. The conclusive 
size of the transferred lepton asymmetry A^ can easily 
be found from the fact that total lepton number per co- 
moving vohimc^ is conserved and that the x-condensate 
has decayed at late times: 

POO 

A^ = -A^lt^oo = / dt- TA^, (8) 
Jo 

where we additionally used equation (7). 

Finally B- and L- violating electroweak sphaleron pro- 
cesses transform parts of the lepton asymmetry in the 



fermionic sector, such that the baryon asymmetry is of 
the same order of magnitude as the lepton asymmetry in 
the DE sector [38]. 



III. QUANTITATIVE TREATMENT 

In order to show that the final asymmetry, which is 
produced in our scenario, can quantitatively explain the 
observed value of the baryon to entropy ratio in the uni- 
verse, it is desirable to get an analytic estimate of the 
asymptotic value Aoo = —A^\t^oo, which we will do in 
the first part of this section. We will then use this esti- 
mate to quantitatively determine the final baryon asym- 
metry of the scenario while we will also present numerical 
results. 

An analytic estimate is indeed possible for a wide range 
of parameters and initial conditions by solving the equa- 
tion of motion (6) of the X-^eld approximately using a 
WKB approach, and then calculating A^ through the 
integral in Eq. (8). Here, we restrict ourselves to the 
parameter space for which the relation 

W,^oc)-.(* = 0)|<^°°d,5JM_«i p) 

is valid. This will be true for the biggest part of the 
considered values due to the growing VEV of the quint- 
essence field. In this case, we can approximate 9 to be 
constant and therefore take the quintessence field to be 
along the real axis without loss of generality.^ Now, the 
equation of motion (6) of the x-ficld can be split into two 
independent linear equations for the real and imaginary 
parts Xi and X2 which are of the form of a harmonic oscil- 
lator with time-dependent damping "/(t) = SH{t) + T(t) 
and frequencies Wi_2(0 = ('^i ± ^2)<j)^{t) + M^, respec- 
tively. For simplicity, we will restrict ourselves to the 
case where 7 <C uji^2 keeping the oscillator from being 
overdamped and assume that the VEV of the quintes- 
sence field gives the main contribution to the mass of the 
X field, fi^ <C (Ai±A2)^o- The requirement for the vaHd- 
ity of the WKB approximation is that the rate of change 
of the frequency and the damping term are much smaller 
than the frequency itself, which in our case simply means 

\m « i-^i- (10) 

For the moment we will assume this equation to hold, 
while the validity of this assumption will be discussed 



^ Let us emphasize that this a very good approximation for the 
dynamics of the x-field, but does not mean that the phase veloc- 
ity could be neglected as far as the dynamics of the whole system 
arc concerned, especially when calculating via Eq. (3) numer- 
ically. The exact parameter space where Eq. (9) and hence the 
WKB approximation is valid is discussed later in this section, 
where we also consider further parameter ranges numerically. 
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Figure 1: Absolute value |x(i)i (solid), and the WKB approx- 
imation (dotted) for Hi„f = (j>o = Xo ~ 10^^ GeV, g = 1, 
Ai = 1, A2 = 0.1, Vo = IO^^pE)"', ki =7 and different val- 
ues of oto ■ For reasons of clearness the plot has been rescaled 
by factors 10~^ and 10~^ for ao — vr/S and 0, respectively. 
For ao 7^ there is a beat frequency u)-/2 = {lji — uji)l2 
corresponding to the large oscillations in \~x\, while the fast 
oscillations of x itself are set by uj+/2 = (cji -\-lj2)/'2. 



later in this section. Now, the result of the WKB ap- 
proximation to leading order in (j>/(j>^ and ^juJx^i is given 
by (see also Fig. 1) 



Xl.2W_ ^1,2(0) Ifeq) 



X?,2 



wi,2(i) 



COS 



dia;i,2(i) 



(11) 

where the initial condition x(0) = has been imposed 
and x(0) — Xoe~'"° is satisfied by the choices = 
Xocos(ao) and = -xosin(Q;o). 

Thus Xi.2(i) are both oscillating, but with different 
frequencies cji.2 if A2 ^ 0. This leads to a modulation in 
the phase velocity of the complex field x and thus of the 
comoving asymmetry, which can be calculated using the 
corresponding expression to Eq. (3) for x: 



^xW = Jsin(2ao)Xo 



/wi(0)w2(0) 

Wl(t)w2(t) 



X cj I sm 



dfw+(?) 



(12) 

with uj± = (jji ± (jj2 — (a/Ai + A2 ± V^i — X2)<p{t). 
Note that only the second part of the damping term^ 
7(t) = 3H{t) + r(t) in Eq. (11), which comes from the 
decay of Xi appears in the asymmetry, while the Hub- 



^ It is now easy to convince oneself that even if different decay 
widths for x ™ different directions in the complex plane were 
assumed the only change in the resulting asymmetry would 
be to replace F in Eq. (12) by the mean value of the corresponding 
decay widths, without changing the main result. 
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Figure 2: The comoving asymmetries A^^^ = n,f,^^{a/ aiif in 
units of [3.2 (pg"')^/^] for (j> (solid), x (dashed) and the WKB 
approximation (dotted, only the leading term proportional to 
UJ+ in the large bracket in Eq. (12) has been kept) for qq = 
— 7r/4 and the other parameters as in Fig. 1. In the two upper 
graphs the x field decays (g = 1, 0.5), the lower graph shows 
the behavior without decay for comparison. In these units, the 
B—L asymmetry nj s is the asymptotic value of A^ for D = 1 
according to Eq. (18). 



ble term gives a contribution (ao/a(t))^/^ which is ex- 
actly canceled. Thus, if there was no decay, F = 0, the 
amplitude of the oscillations in the (comoving) asymme- 
try would stay constant, whereas the frequency increases 
since the VEV of the quintessence field </> grows when it 
rolls down its potential, see the lower graph in Fig. 2. 
When F > 0, the increase in the frequency still occurs, 
but the amplitude is damped by the exponential term in 
Eq. (12), see the two upper graphs in Fig. 2. Further- 
more, the asymmetry is proportional to sin(2Q;o), i.e. it 
is zero if there is no initial phase difference (ao = 0), and 
its modulus is maximal if — ±7r/4. 

It is instructive to consider the case when A2 <C Ai , be- 
cause then uj\ K. 4Ai0^ ':$> ujt ~ (A|/Ai)0^ and the first 
term in the large brackets of Eq. (12) dominates, giv- 
ing the leading oscillations of with a "low" frequency 
cij- proportional to the parameter A2, superimposed by 
a "fast" oscillating component whose amplitude is sup- 
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pressed by a relative factor lo-/uj+ (see Fig. 2). 

Finally, the analytic expression (12) together with a 



decay rate of the form V{t) ~ g j {%'K)\/\\(^{t) can be 
used to calculate the total asymptotic asymmetry (per 
comoving volume) via the integral in Eq. (8): 



^oo = 2 sin(2ao) Xq 00 



Stt r2. 



- L>J_ 



1 - 



(13) 

Note that the time-dependence of 0(i) disappears in the 
final result and only the initial values 0o and xo en- 
ter as dimensionfuU quantities. Furthermore, A^^ van- 
ishes if either F cx g^/iSiT) or A2 oc uj+lo^ 
since in these limiting cases no transfer of asymmetry 
between x a-nd the SM, or, respectively, x and the quint- 
essence sector would be possible. If, on the other hand, 
(/2/(87r) « A2/A1 <C 1, both transfers are of compa- 
rable "strength" and the asymmetry is simply A^o ~ 
sin(2ao) Xq •\/^0o/2, independent of 5 and A2. For given 
initial values (pQ, Xoi Q^o and fixed Ai, this is the maximum 
asymmetry that can be transferred. 

It turns out that Eq. (13) gives a very accurate estimate 
for in the case 



(14) 



since conditions (9) and (10) are both fulfilled from the 
beginning (i = 0). This in return can easily be seen from 
the approximate analytic solution for the quintessence 
field for early times (when 4){t) ^ Mpi) 



(15) 



In this domain the asymmetry is independent of Vq , but 
(for fixed xo) only depends on the initial quintessence 
VEV (jiQ (corresponding to the vertical contour lines in 
Fig. 3). 

If, on the other hand. 



\v'{M\:>4>l xl 



(16) 



condition (9) on the phase is actually also satisfied, but 
condition (10) is only fulfilled for times t > ^match = 
(4/|y (^o)!)^'''^ as can again be easily seen using (15). To 
get an estimate of the asymmetry in this case one can use 
the early-time solution (15) for t < imatch and then match 
the WKB solution for t > tmatch- Since (^oimatch ^ 1 (see 
Eq. (16)) the x-field can be approximated to stay static 
until tmatch, and thus the only modification in the WKB 
results obtained above, especially in the final asymmetry 
(13), is to replace the initial time t — Ohy t — imatch and 



00 -> 0(imatch) ~ VWiMl = ^2aK)/Mpi. (17) 

In this domain, in contrast to the previous case, the fi- 
nal asymmetry is thus independent of the initial quintes- 
sence VEV 00 and only depends on the initial potential 
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Figure 3: This contour plot shows the final asymmetry factor 
n from Eq. (18) for the same parameters as in Fig. 2 (g = 1) 
apart from varying values of the initial quintessence VEV (j>o 
and the initial potential energy Vo- The light dot-dashed lines 
divide the plot into three domains: In the right part, due to 
Eq. (14), 00 ^ Xo0O) l^'(0o)|, the asymmetry depends on 
(j>o but not on Vo, whereas well inside the upper left part it 
is vice versa since \V'{(f>o)\ 2> <l>Oy Xo from Eq. (16). This 
is in accordance with the WKB results (20) inside each of 
these domains. In the remaining part in the lower left comer 
backreaction effects of x on (p lead to more complex dynam- 
ics. Note that n/s = kD with D — 1 . . . 10~®, see Eqs. (18) 
and (19). 



energy'* Vq (corresponding to the horizontal contour lines 
in Fig. 3). 

To determine the final baryon number, we need to 
calculate the ratio of the produced B-~L asymmetry 
n = —Aoo{ao/a)^ and the entropy density [39], 



^9 iTK)T^[-) ~3.2pj, ^-j. 



where Tr, = g* (T[{)tt'^T^/30 and apj are the temper- 
ature, energy density and the scale factor at reheating, 
respectively, and the number g*(TR) of active degrees of 
freedom has been assumed to be of order one hundred. In 
principle, the final asymmetry of the scenario is indepen- 
dent of the physics between inflation and reheating, how- 
ever, inflationary dynamics can enter indirectly through 
the expansion rate when relating the corresponding scale 
factors and critical energy densities via 



where u!{t) is the equation of state parameter of the dom- 
inant component, i.e. the inflaton and later radiation. 
Thus the final asymmetry relative to the entropy density 



For 4>o <C Mpi we have V{4)o) = Vo. 
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is given by^ 



3.2(p§"*)3/4 



exp 



dt{l - 3u{t))H{t) 



iD 



(18) 

The first factor k depends only on the dynamics of the 
quintessence and the mediating field shortly after infla- 
tion, whereas the second factor D describes the dilution 
of the asymmetry between the end of inflation and the 
time when w(i) ^ 1/3. In fact, l/D is equal to the 
spatial expansion of the universe after inflation relative 
to a radiation dominated cosmos. Note that the upper 
limit in the integration can be formally extended beyond 
reheating, since from then on w = 1/3. This makes it 
obvious that any inflationary scenario can only have an 
impact on the final baryon number through its equation 
of state and not through further characteristics, e.g. the 
reheating temperature^. 

To get a plausible range for D we will consider various 
possibilities: for example, if the inflaton field / oscillates 
in a potential of the form /™ with even m, its mean equa- 
tion of state (averaged over a whole oscillation period) is 
Lo — — const [40] . Assuming that this is true for the 
whole time between the end of inflation and reheating we 
get D = (2.4rR/(pg"*)i/4)(4/'"-i). For Tr = lO'' GeV, 
-ffinf =10^^ GeV and m = 2 this gives D ^ 10^^', whereas 
for m = A there is no dilution at all yielding D = 1, with 
or without preheating [40]. Another example is a model 
with preheating and a quadratic minimum [m = 2) for 
the inflaton. For the case recently discussed in Ref. [41], 
one finds D w 10""'^ . . . lO"''. Here we will adopt the point 
of view that the dilution factor D will be roughly in the 
range 1 . . . 10~^, yielding a range for k 



10-1" < K < 10" 



(19) 



that is compatible with the observed value 
n/s ~10-i° [1]. Using the analytic estimates (13) 
and (17) one obtains 



K w — -— sin(2ao) 



Xo 



<t>a 



7.3H,, 



Mp, 
1/3 



3/2 



1 (s.iVoY'-' f m„! 



7/6 



(20) 



withj\r= ^-^^^^il-uj^_/ujl). The upper and lower 

line of Eq. (20) correspond to the cases (14) and (16), re- 
spectively. For coupling constants g^/{S'K) ~ A2/A1 -C 



^ Since only a part of this asymmetry is transferred from the lep- 
tonic to the baryonic sector tiicrc will be an additional sphaleron 
factor of order one, see Rcf. [38] . 

^ Of course, there arc inflationary scenarios in which the time of 
reheating also marks a change in the equation of state, in which 
case the amount of the produced baryon asymmetry and the 
reheating temperature are not completely independent. 



1 ~ Ai, A/" is of order one. Generically assuming 
sin(2ao) ~ -1 and xo ~ -H"inf ~ 10^^ GeV - 4 • lO^^Mpi 
we thus find values for k of roughly lO^^^^o/^inf and 
10-^(Vb//Oo'^'*)i/^ respectively, which lie inside the al- 
lowed range (19) for a huge variety of initial conditions 
of the quintessence VEV and its potential energy after 
inflation (see Fig. 3). 



IV. FERMIONIC COUPLINGS 

There are various possibilities to couple the mediating 
field X of our toy model to the fermionic sector for which 
we will give two examples in this section. 

Since the whole scenario is B—L conserving, the right- 
handed neutrinos do not have a Majorana mass term and 
carry lepton number one. As the x-field is also a gauge 
singlet it can be coupled to one (or more) right-handed 
neutrino (s) A'' via 



gN^^Nx + h.c. 



(21) 



This would assign the desired lepton number to the x 
field and implies the given decay term Fx in the equa- 
tion of motion (6). However, if this was the only addi- 
tional coupling, the asymmetry would remain stored in 
the right-handed neutrino sector due to their tiny Yukawa 
couplings [6] . In this case one would have to find a way to 
transfer the lepton asymmetry to the SM sector, where 
sphalerons could create the final baryon asymmetry. This 
could be realized by a light scalar field S (with mass 
-f^inf ^ "'-s ^ TeV) with the same quantum numbers 
as the SM Higgs field and a quadratic potential. If this 
scalar would start with a typical VEV of order i?i„f at the 
end of inflation, it would remain frozen in until H drops 
below the mass of the scalar. In this case an additional 
coupling term 



g' mS + h.c. 



(22) 



with a left-handed lepton doublet £ and suitable cou- 
pling constant g' should lead to the required left-right- 
equilibration of the neutrinos as H becomes smaller.^ 
When H drops below the mass of the scalar the VEV 
of S would start oscillating and finally vanish. It might 
even be possible to avoid the introduction of such a new 
particle by introducing a similar coupling between the 
inflaton and the neutrinos. In this latter case one would 
get additional constraints for the inflationary scenario. 

One possibility, which works without the introduction 
of additional particles while also being independent of 
the inflationary model, would be to couple the x field 
directly to the left-handed neutrinos v via the term 



g v'^ux + h.c. 



(23) 



^ The quartic potential of the SM Ifiggs and its small Yukawa 
couplings to neutrinos in the case of Dirac neutrinos seem to 
prevent it from being a suitable candidate for this job. 
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Due to the heavy mass of x this should be consistent with 
observations. However, the exphcit breaking of gauge 
invariance makes this option less attractive. 



V. STABILITY AND Q-BALLS 

In this section we discuss possible sources of instability 
in our model. For instance, particle processes and spatial 
inhomogeneities could disturb the scalar field evolution 
described in the previous sections. One peculiarity of 
complex scalar fields with a conserved global quantum 
number is the potential formation of so-called Q-Balls 
[42] . A Q-Ball is a non-topological soliton and represents 
a spatially extended scalar field configuration carrying a 
non-zero charge associated to the global symmetry. For 
our scenario it is essential that the complex quintessence 
field does not promote the formation of Q-Balls, which 
could make the field unstable and inhomogeneous and 
thereby destroy its DE properties. To analyze the stabil- 
ity of the scalar fields in our model one has to to solve 
the corresponding equations of motion including the spa- 
tial derivatives that were omitted in Eqs. (5) and (6). 
Since the numerical solution of this system of partial dif- 
ferential equations is beyond the scope of this paper (see 
e.g. [43]), let us instead proceed with a much simpler ap- 
proach that might give us at least some hints towards the 
stability properties of our model. Since we also only want 
to show that instabilities should not be a problem for our 
model in general, we only consider the initial conditions 
and parameter values as in Fig. 1 {g = 1) and neglect 
perturbations in the metric. 

Following Refs. [44, 45] we study the evolution of 
small inhomogeneous perturbations in the scalar fields. 

Substituting (b{x,f) 0(t) + 5(j){x,f) and xi^^'^) ~^ 
x{t) + 6x{x, t), where the first parts denote the homoge- 
neous solutions of the previous sections, we can derive the 
equations of motion for 6(j) and 6x from the Lagrangian 
( 1 ) . Keeping only the terms linear in the small perturba- 
tions, we find a mass mixing matrix for S(f> and 6x, where 
the dominating contributions on the main diagonal, at 
early times, are of the order y(|(/)|)/|^| and Ai|0p for 
S4> and Sx, respectively, while the off-diagonal terms are 
of the order AiX(^. As we checked numerically, the off- 
diagonal entries never become solely dominant, therefore 
we will consider the perturbations of (p and x separately 
in this simplified analysis. 

Now, we write the complex scalars (f> and x the 
form (p — ¥'6*^, X — lx|c*'^ with real functions (p, 
and a. Within the quintessence field (p we consider per- 
turbations of the type 6(f = (5v?o exp(r2(t) -|- ikx), 59 = 
d9o exp{Q,{t)+ikx), that are characterized by a wavenum- 
ber fc = |fc| and a time-dependent function Cl(t). From 
the Lagrangian (1) we derive the equations of motion for 



the quintessence field 

= <f + 3Hip-a-'^V'^ip-ipe'^ + a-'^ip{S/9f 
+ V''((^) + |xp(^[Ai + A2Cos2a], (24) 

= e + 3He~a-^v^e + 2e^ 

- 2a-'^{\7ip){V9) - X2\xfsm2a, (25) 

where a = 9 — a is the phase difference of both scalar 
fields. By substituting (f — > ip + Sip, ^ — > 6 + 69 and keep- 
ing only terms at most linear in the perturbations 5(p, 59, 
one obtains a system of equations for S(p and 59. The 
condition for this system to have a non-trivial solution 
reads 

= (n + n^ + 3Hn + ^-9'^ + v"{ip) 

+ lxl^[Ai -FA2Cos2a]) x (O-^fi^ 

+ 3Hn + ^ + 2(l^ - 2X2\xf cos2aj 

+ 49 (n9 + X2\x\^ sm2a^ x (fi-^). (26) 

If this equation for ri(t) has any solution with Re(f7) > 
for a period of time, then the perturbation mode 
~ exp(f2(t) + ikx) is growing exponentially, possibly in- 
dicating an instability in the scalar field (p, which could 
lead to the formation of Q-Balls. The following discus- 
sion is divided into three parts. First, during the early 
moments of evolution, when the mediating field x has 
not yet decayed, we will study the stability of the model 
mainly numerically. The second and third parts concern 
the time directly after the x-decay and the very late cos- 
mological era, respectively, where we can now treat the 
problem completely analytically. 

Starting at the end of inflation we solve Eq. (26) 
numerically for the function Cl{t) and several different 
choices of the wavenumber {k/a)"^ and the initial phase 
difference ao- In some cases one can observe positive 
values of Re(f2), which appear maximal for small k. As 
an example the fc = case is shown in Fig. 4. How- 
ever, for these unstable modes, Re(f2) turns out to be 
much smaller than the Hubble scale ffinf so that there is 
not enough time for the instabilities to develop. In ad- 
dition, a wavenumber smaller than -ffinf implies that the 
wavelength of the perturbation is larger than the Hubble 
radius and its evolution is therefore suppressed by the 
cosmological expansion. For modes with large wavenum- 
bers (fc/a)^ S> -ffj^^j the numerical solutions of Eq. (26) 
do not show growing perturbations, as expected. 

To study the late-time behavior of the inhomogeneities, 
we simplify Eq. (26) by neglecting <C i.e. as- 

suming a slowly changing uj = tl since the background 
solution is now slowly varying. In addition the mediating 
field X has already decayed (]x| = 0). Thus, we obtain 

= (w'^+2.Hw+^-9'^ + V"{ip)) 

X (a;2^3ifa;+f^ + 2u;|) +4a;^2(a;-|),(27) 
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Figure 4: The real part of Q{t), where the function n{t) 
describes the evolution of the inhomogeneous perturbations 
Sifi, 56 oc exp(f2(t) -f ikx) for different values of the initial 
phase difference: ao = 7i"/4, tt/S, 7r/16. For this numerical 
solution the wavenumber (k/a) and the initial value of Q. are 
taken to be 0. The remaining parameters are the same as in 
Fig. 1. 



where solutions with a positive real part of uj might in- 
dicate an instability. We find that the range of val- 
ues for the wavenumber k/a that solves this equation 
for Re(a;) > is given approximately by 



< 



- 1 < 

a 



= e''-v"{>p). 



(28) 



Hence, the late-time evolution of the quintessence 
field is stable for 9'^ < V"{ip). To verify this 
condition let us assume for the comoving asym- 
metry of the quintessence field a maximal value 
of \A^\ = 2\e\ip^{alaQf - 10-«(p^"*)3/4, where the ini- 
tial critical energy density is given by p"'' ~ H'^^^M^y 
Although Lp starts with the value i^jq ~ ^^inf, it quickly 
reaches the Planck scale Mpi, therefore we take the mod- 
erate value (/? ' 



Vi/infMpi at early times so that 



1^1 



/ crit\3/4 
-in-S \Pq ) 
9 



< 10-V^infMpi. (29) 



Even if we now choose the low value Vq ~ Hf^fMpi in the 
quintessence potential V{ip) = Voe~^^'^^-'^^\ we obtain 



< 10- 



IQ-^ < 1, 



(30) 



where we used iJinf ~ 10 ^ Mpi and ^ 10. At very late 
times in the cosmological evolution with a low Hubble 
scale H < i?i„f we find V"{ip) = 0{H'^) [46], ip > Mpi 
and {a/ao)~^ — (H/Hi^f)^/^ for a long period of radia- 
tion dominance. Thus, the quintessence field is also sta- 
ble in this epoch: 



V"{ip) 



10 



-16 



H 



Ah 



< 1. 



PI 



(31) 




X 



Figure 5: Example for an additional process between (p and 
the leptonic sector. 



We conclude that the quintessence scalar field in our sce- 
nario does not seem to suffer from instabilities, at least 
within the context of this simplified analysis. In prin- 
ciple, we could perform a similar analysis also for the 
mediating field x- However, since it decays into neutri- 
nos very quickly, a possible instability of x should not 
harm our scenario. Additionally, Q-Balls built from this 
field could also decay into leptons [47-49]. For a deeper 
discussion of Q-Balls see e.g. Ref. [50] and the references 
therein. 

In the last part of this section we want to motivate the 
statement that particle processes only enter the equa- 
tions of motion of the fields in form of the decay term Tx 
in Eq. (6) and make it plausible that this might even 
hold within a full quantum-mechanical treatment of the 
quintessence field. However, one has to admit, that such 
a treatment of the quintessence field would go far beyond 
the scope of this paper and more subtleties would arise 
than the ones mentioned here. 

Still, we notice that one could also allow additional par- 
ticle processes like the one in Fig. 5, which might desta- 
bilize the (f) condensate. However, even if this was done 
it is unlikely that they would have a large impact on the 
scenario due to the following arguments. We first notice 
that at the beginning of the scenario the masses of the 
scalar particles (proportional to the second derivative of 
the corresponding potential) are of the order of Hi^f or 
smaller. By dimensional analysis the reaction rates Fi 
for the mentioned processes should therefore be of the 
same order of magnitude at this time (assuming coupling 
constants to be of order one). However, before the cor- 
responding time elapses the mass of % has increased by 
several orders of magnitude. Since these masses appear in 
the denominator of the reaction rates the ratio Ti/ H will 
now always be smaller than one and the reaction should 
never have enough time to take place as also resonances 
do not seem to occur. 

We also note that additional decay channels of x would 
not harm our scenario as long as the lepton asymme- 
try is transferred to SM particles before the electroweak 
phase transition. A decay channel which could poten- 
tially be dangerous is a decay into hypothetical quanta 
of the quintessence field, if they represent physical de- 
grees of freedom, via the couplings of the form \i^24'^X^- 
However, a simple estimate for the decay rate T2x-^2(t> ~ 



9 



•^1 2lxP/(87i'™x) shows that it is both suppressed by 
the increase of and the decrease of |x| as compared 
to the rate ^x^'^i' — Q^'^xl l^"") fo'" decay into neu- 
trinos. 



VI. SUMMARY AND CONCLUSIONS 

In this work we have proposed a B—L conserving 
baryogenesis scenario, in which a B—L asymmetry is hid- 
den in the DE sector. This is achieved by the introduc- 
tion of two complex scalar fields carrying Icpton number, 
where one of them is responsible for DE in a quintessence- 
like manner while the other one mediates between this 
field and the SM particles (including right-handed neu- 
trinos). An initial phase difference between the VEVs of 
these two fields results in compensating but non-zero Icp- 
ton asymmetries for them. While the quintessence field 
has a tracking behavior, the mediating field acquires a 
large mass and eventually decays, hereby transferring its 
asymmetry to the fermionic sector. Finally sphalerons 
partially transform this lepton asymmetry to a baryon 
asymmetry. We quantitatively determined the resulting 
baryon asymmetry numerically and also gave an analytic 
formula based on a WKB approximation. Both results 
are consistent and show that the observed value can eas- 
ily be achieved. Furthermore, wo have analyzed potential 
sources of instability like Q-Ball formation and did not 
find any hints for an unstable behavior. 

Of course, a deeper understanding of several topics is 
desirable. A thorough knowledge about the VEVs of 



scalar fields after inflation, especially in the case of com- 
plex potentials, could help to restrict the parameter space 
of our model immensely. Also, as stated before, a full 
quantum-mechanical treatment of the quintessence field 
is wish-full. This could enable us to further restrict the 
couplings in our scenario, while it might also make way 
for new baryogenesis models in which the quintessence 
field is directly coupled to SM particles. 

The possible B—L quantum numbers of the vacuum 
could also have interesting consequences for additional 
couplings between DE and SM fields or dark matter. 
Moreover quantum numbers of the quintessence field 
could include hints at its origin. 

In conclusion, we presented a specific model that leads 
to the observed baryon asymmetry of the universe with 
a strict conservation of B—L and the storage of a corre- 
sponding asymmetry in the DE sector. 
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